arXiv:1507.03365vl [cond-mat.mes-hall] 13 Jul 2015 


Direct surface charging and alkali-metal doping for tuning the interlayer magnetic 
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The continuous reduction of magnetic units to ultra small length scales inspires efforts to look for 
a suitable means of controlling magnetic states. In this study we show two surface charge alteration 
techniques for tuning the interlayer exchange coupling (lEC) of ferromagnetic layers separated by 
paramagnetic spacers. Our ab — initio study reveals that already a modest amount of extra charge 
can switch the mutual alignment of the magnetization from anti-ferromagnetic to ferromagnetic 
or vice verse. We also propose adsorption of alkali metals as an alternative way of varying the 
electronic and chemical properties of magnetic surfaces. Clear evidence is found that the interlayer 
magnetic order can be reversed by adsorbing alkali metals on the magnetic layer. Moreover, alkali 
metal overlayers strongly enhance the perpendicular magnetic anisotropy in FePt thin films. These 
findings combined with atomistic spin model calculations suggest that electronic or ionic way of 
surface charging can have a crucial role for magnetic hardening and spin state control. 


I. INTRODUCTION 

Efficient ways of manipulating memory units and logic 
devices became an emergent motive in the field of spin 
electronics^’^, in addition to minimizing their sizes. In 
the last few decades the reading scheme of the modern 
hard disc drive was reduced to the nanoscale level fol¬ 
lowing the discovery of giant magnetoresistance effect in 
metallic multilayers by Grunberg and Fert^’^. Current 
technology uses tunneling magnetoresistance for reading 
the state of magnetic units^’®. It relies on the relative 
alignment of the magnetic layers interspaced by coupling 
medium, which results in the difference of the tunneling 
resistance. The existence of such intriguing characteris¬ 
tics is related to the delocalized electrons within the cou¬ 
pling medium that are confined. Apart from the advance¬ 
ment in the reading mechanism of the state of the mag¬ 
netic units, the writing head still makes use of magnetic 
field which has less capability of controlling magnetic 
states on short length scales. Alternatively, such mag¬ 
netic states can be reversed using the current-induced- 
spin transfer torque technique in which the spin-polarized 
current tunnel through the junction and consequently 
producing spin-torque^’®. It has also been shown that 
the mutual magnetization in magnetic tunneling junc¬ 
tion can be tuned by applying a voltage and this effect 
has similar order of magnitude like that of spin transfer 
torque®^^^. 

Recently a more efficient method has been revealed, 
on the basis of magnetoelectric coupling, for tuning the 
spin alignment in a ferromagnet (FM) by both Density 
Functional theory^®^^® and experimental studies^’^^’^®. 
For instance, the anisotropy properties of 2D ferromag¬ 
net could be tuned either by varying its intrinsic charge 
carriers^^’^® or the oxidation state by electric means^®. 
On the other, hand external electric field was used in 
combination with multiferroic materials to reverse the 
relative alignment of ferromagnets separated by non¬ 
magnetic medium (FM/NM/FM)®®“^^. Theoretical stud¬ 
ies have shown that the interlayer exchange coupling and 


magnetoresistance effect are inter-related^®’^^. Nonethe¬ 
less, direct manipulation of the interlayer exchange cou¬ 
pling is one of the feasible and perhaps the least investi¬ 
gated approaches for a tuning magnetic order. It would 
even be more compelling to find a possible way of tuning 
magnetic coupling or anisotropy in planar nanostructures 
with ionic doping. 

Our study verifies a possible route for manipulating 
the lEC and the relative alignment of the magnetic lay¬ 
ers interspaced by non-magnetic multilayers, employing 
direct surface charging and alkali metal adsorption. The 
former mechanism relies on sequential variation of the in¬ 
trinsic charge concentration of FM/NM/FM multilayer, 
and the results related to this phenomena are presented 
in Sec. III. Additionally, in Sec. IV we propose that the 
process of surface charge alteration can be accomplished 
by adsorbing free-electron systems, alkali metals, on sur¬ 
faces. Indeed, this mechanism alters the electrochemical 
features of magnetic surface and leads to reversal of the 
relative spin alignment. Apart from tuning the inter¬ 
layer magnetic order the magnetic anisotropy of Fe/Pt 
multilayer can be significantly enhanced by using alkali 
metal deposition. To better appreciate such increase of 
the MAE, in Sec. V the hysteresis loop for FePt mul¬ 
tilayer is plotted through the atomistic spin model ap¬ 
proach. The discussion and analysis on these findings is 
presented in Sec. VI, followed by conclusion in Sec. VII. 


II. COMPUTATIONAL DETAILS 

A thorough ab —initio study was performed within the 
projector augmented wave technique (PAW)®® as imple¬ 
mented in the VASP code®®’®^. The local spin density 
approximation (LSDA) is used for the exchange and cor¬ 
relation interactions®®. Plane wave basis set is used to 
describe the Kohn-Sham wave function, and in all calcu¬ 
lations the plane wave cut off energy of 430 eV was used. 
A dense /c-point mesh of 21x21x1 and smearing width of 
2 meV, that are tested for their reliability, were employed 
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to obtain a more accurate value of MAE and lEC. We 
have used a supercell approach in which two layers of Fe 
separated with finite number of Cu spacers are supported 
by Pt(OOl) substrate^®. The surface charging effect is in¬ 
troduced by sequentially varying the number of valence 
electrons. Additionally, a uniform charge background is 
assumed^®. Doping of alkali metals is done by capping 
the top magnetic layer with two different structures of 
Na or Li overlayers, namely p(lxl) and c(2x2). For the 
latter configuration a 2x2 supercell is used. The pairwise 
exchange coupling or simply the interlayer exchange cou¬ 
pling (Eiec) was calculated as Ejec = Eaf — Efm, and 
for all charge states the spins were aligned along their 
easy axes. For 2x2 supercell the interlayer coupling is 
calculated as, Eiec = j[Eaf — Efm]- To determine the 
easy axis of magnetization we have performed a fully- 
relativistic calculation including spin-orbit coupling^^. 
The magnetic anisotropy energy (MAE) is evaluated by 
taking the energy difference between two axis of magne¬ 
tization, parallel ([ 100 ]) and perpendicular ([ 001 ]) to the 
surface plane. Similar procedure is followed to calculate 
the MAE for antiferromagnetic configuration and the two 
spins are oriented anti-parallel to each other for both axis 
of magnetization. Here after, we shall denote the charge 
added to the multilayers (net charge) as (?+(-) whereas n 
stands for the charge density of the entire system. 

In order to study the magnetization reversal in FePt 
planar nanostructures, we use atomistic spin model 
approach.It is essential enough to investigate this phe¬ 
nomena by mapping the magnetic properties already in¬ 
vestigated on the basis of electronic structure theory with 
atomistic spin models. In the latter case the energetics 
of interacting spin with an effective magnetic moment is 
described by spin Hamiltonian which has a form 


H = •S, -fc„^(S,-e)® 

i 

where the first, middle and last terms describe the ex¬ 
change interaction {H^xc), magnetic anisotropy (Hmae) 
and the external magnetic field (Happ) or Zeeman term 
respectively. Jij is the exchange coupling parameters be¬ 
tween atomic sites i and j, and the local spin moments 
are denoted with unit vector Si and obtained from 
the realistic atomic moments as Si = ms/\ms\. The ex¬ 
change coupling parameters and the single ion anisotropy, 
in Eq. 1, are obtained from calculations on the basis 
of density functional theory. The magnetization curves 
are simulated by solving the stochastic Landau-Lifshitz- 
Gilbert (LLG) equation. At the atomic level the LLG 
equation is written as 


dt 




where 7 = 1.76 x 10^^ T“^s“^ is the absolute value 
of the gyromagnetic ratio, A is the damping parameter 
and the atomistic effective magnetic field is denoted with 
The effective magnetic field is derived from the 
spin Hamiltonian shown in Eq. 1.^^ The FePt multilayers 
are modeled with two Fe layers separated with an effec¬ 
tive medium and it extends as 6 nm x 7 nm within the xy 
plane. When simulating the magnetization curve thermal 
effect and magnetostatic field are taken into considera¬ 
tion that are included in the effective magnetic field^^’^^. 
Further discussion on the hysteresis loop of FePt multi¬ 
layers is presented in Sec. V. 


III. SWITCHING INTERLAYER EXCHANGE 
GOUPLING WITH SURFACE CHARGING 


First, we present the results concerned with the influ¬ 
ence of direct surface charging on the relative magneti¬ 
zation of two magnetic layers separated with Cu spacer. 
When the charge neutrality of metallic multilayers is per¬ 
turbed by additional positive or negative charges, it leads 
to spatial charge re-distribution mainly towards the sur¬ 
face. Whereas, a reduced variation of the spatial charge 
is also observed underneath the surface. Such charge re- 
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FIG. 1. (Color online)(a) The interlayer exchange coupling 
(lEC) between the magnetic layers in Fe/CuN/Fe multilayers 
on Pt(OOl) as function of the net charge in the system. As an 
illustration such relations are plotted for Cu spacer thickness 
of N = 4 (circle) and N = 7 (triangle), (b) A schematic 
description of Fe/Cu 4 /Fe and the relative alignment of the 
magnetization as result of charge injection or removal. 

distribution can even be spin-dependent which can alter 
the magnetic order^^’^®. Meanwhile, the exchange in¬ 
teraction across the coupling medium in Fe/GuN/Fe tri¬ 
layer, supported on a Pt(OOl) substrate, determines the 
magnetic state to be either ferromagnetic (EM) or anti¬ 
ferromagnetic (AF). For such multilayers we demonstrate 
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a strong effect of direct surface charging on the relative 
magnetic order. As an illustration, in Fig. 1 the variation 
of the interlayer exchange coupling for Fe/CuN/Fe multi¬ 
layers in response to the addition or removal of charges is 
depicted for Cu spacer thickness of N = 4, and N = 7 lay¬ 
ers. The Fe layers separated with four layers of Cu (un¬ 
charged) are coupled antiferomagnetically, and the Eiec 
is found to be -2.5 meV. The antiferromagnetic coupling 
can be further enhanced by negative charge doping. As 
an example, by injecting negative charge of the order of 
Aq = 1.2 e per unit cell^^ Fe/Cu 4 /Fe, the absolute value 
of the Eiec is increased to 12 meV. The most fascinat¬ 
ing finding is observed when the electron concentration 
is reduced from the neutral Fe/Cu 4 /Fe, which results in 
the switching of the relative magnetic order from anti¬ 
ferromagnetic coupling to ferromagnetic one. Explicitly, 
switching from AF configuration to FM takes place when 
the system is charged with 0.2 h per unit cell, and further 
increment of the extra holes to 1.2 h per unit cell leads 
to Eiec of ^ 5 meV. 

An essential issue now could be the interplay of finite 
size effect, i.e. thickness of the Cu spacer, surface charg¬ 
ing and the interlayer exchange coupling. In order to in¬ 
vestigate this phenomena the exchange coupling (Eiec) 
with respect to charge injection is investigated for thicker 
Cu layers. In Fig. 1 (a) we have employed seven layers of 
Cu as mediator in which the lEC of the uncharged system 
is 3.6 meV. Similar to the case of four layers, here space 
charge variation has significant effect on the lEC. Addi¬ 
tion of certain amount of negative charges switches the 
magnetic coupling to AF. Interestingly, for Fe/Cuy/Fe 
mutilayers a transition of the easy axis of magnetization 
from in-plane to out of plane directions and the reversal 
of the relative magnetic order are observed at the same 
charge state, i.e. at 0.6 e/unit cell. In another multilayer, 
i.e Fe/Cug/Fe trilayer, we found that the magnetic mo¬ 
ments are counter-aligned for all charge states. But, the 
influence of negative charge injection on the lEC is much 
stronger than the positive ones. Precisely, the absolute 
value of the lEC for the neutral system (^6 meV) is re¬ 
duced by 75 %, as soon as the system is injected with 
1 e per unit cell. One should note that for the experi¬ 
mental realization of effective surface-interface charging 
the alloy multilayers should be electrically decoupled, by 
semiconductor or insulator as shown in Fig. 1 (b), and 
hence the net charge is trapped in the upper layers. 


in nanostructure can be achieved by electrolyte^®^^^, or 
by adding dopants^®, whereas using scanning tunneling 
microscope or electron force microscope tip is also in¬ 
dispensable option®®. Doping of magnetic entities with 
alkali metals or halogens can lead to a significant varia¬ 
tion of the electrochemical properties'®’"^® as well as their 
spin properties"®®’'®®. Moreover, alkali metals are illustra¬ 
tive free-electron systems and have relatively small ion¬ 
ization energies. Here we verify that alkali-metals de¬ 
posited on magnetic surface can be employed as a nat¬ 
ural way of varying the electronic and chemical proper¬ 
ties of FM/NM/FM trilayers and the magnetic properties 
could be controlled therewith. Particularly, the interplay 
between the deposition of alkali metals and the change 
in the interlayer magnetic coupling is investigated. In 
Fig. 2 (a) the values of the lEC for Fe/CuN/Fe and 
Na/Fe/CuN/Fe multilayers is presented. The magnetic 
coupling of Fe/CuN/Fe oscillates between FM and AF 
order as the thickness of the spacer increases, which is 
related to the quantum well states (QWS) in the Cu 
multilayers®®’'®®. In fact, the results for six, seven and 
eight Cu layers of the neutral systems are compared 
with available studies based on Korringa-Kohn-Rostoker 
method®®, and the magnetic order as well as the exchange 
coupling energies are in a very good agreement. 



Number of Cu layers (N) 


IV. ALKALI-METAL-INDUCED SWITCHING 
OF THE INTERLAYER MAGNETIC ORDER 

So far we have shown the impact of charge injection on 
magnetic order by controlling the intrinsic charge car¬ 
riers. Albeit, controlling the number of charge carri¬ 
ers in a slab could be experimentally challenging task. 
Thus, the efforts looking for an alternative method of 
varying the surface electronic structure are underway. 
Recently it was shown that incorporating a net charge 


FIG. 2. (Color online) (a) Impact of Na adsorption on the 
interlayer exchange coupling (lEC) in Fe/CuN/Fe (triangles) 
and Fe/CuN/Fe/Na- p(lxl) (squares) for different thickness 
of Cu spacer. The values of the lEC for c(2x2) configura¬ 
tion of Na (top) on the magnetic layer are represented with 
circles. The inset shows the p(lxl) and c(2x2) structures 
of Na or Li adsorbates on the magnetic surface, (b) Sim¬ 
ilar plot on the variation of the lEC for Fe/Cux/Fe and 
Fe/CuN/Fe/Li (rectangles and circles) multilayers with re¬ 
spect to Cu spacer thickness. The lEC of Fe/CuN/Fe/Li 
multilayer with p(lxl)[rectangles] and c(2x2) [circles] config¬ 
uration of the Li capping layer is also presented. 
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We have considered two structures of Na (Li) on the 
Fe surface, namely Na (Li)-p(lxl) and Na (Li)-c(2x2) 
structures, in order to have more insight on the influ¬ 
ence of the Na (Li) coverage on the lEC. The former 
and latter adsorbate structures covers the magnetic sur¬ 
face fully and partially, respectively. The configuration 
of these two structures on the magnetic surface is shown 
in Fig. 2 (a) (see the inset). The Na ions are adsorbed at 
the hollow site which was identified as the most favorable 
one for alkali adsorbates^®’"*®. Regarding the magnetic 
coupling, adsorption of Na ions on the magnetic surface 
brought an intriguing change to the interlayer exchange 
coupling. As shown in Fig. 2 (a) (square), when the Fe 
layer is fully covered with Na [p(lxl)] all the magnetic 
states of Fe/Cux/Fe is switched from FM to AF or from 
AF to FM, except for N = 7. Especially a strong change 
on the lEC, accompanied by magnetic reversal, is ob¬ 
served when the the Cu spacer thicknesses are four and 
eight. Reducing the Na layer coverage by half, from Na- 
p(lxl) to Na-c(2x2), can also lead to significant change 
of the lEC. As shown in Fig. 2 (a) (circles), capping the 
Fe/CuN/Fe multilayer with Na-c(2x2) changes the mag¬ 
netic order from FM (for the uncapped multilayers) to 
AF, when the spacer thicknesses are five and seven Cu 
layers (N = 5 and N = 7). Moreover, the same config¬ 
uration of Na varies the interlayer magnetic coupling in 
Fe/Cu6,8/Fe multilayer from AF to FM one. But, cov¬ 
ering the top Fe layer of Fe/Cu 4 /Fe with sub-monolayer 
Na-c(2x2) does not switch the magnetic order. Switching 
relative magnetic order can be accomplished by increas¬ 
ing the Na coverage, suggesting that the lEC depends on 
the percentage coverage of Na overlayer. Additional man¬ 
ifestation of the effect of Na coverage on the lEC can be 
seen on Fe/Cuy/Fe trilayer. In this multilayer when the 
Na coverage is reduced by 50 %, from p(lxl) to c(2x2), 
the lEC changes from 2 meV (FM order) to -7 meV (AF 
order), respectively. These observations lead us to the 
conclusion that variations of the lEC induced by Na de¬ 
pends on the configuration of the adsorbate (percentage 
coverage) as well as on the spacer thickness. 

Similar investigations have been done by using Li ions 
as dopant instead of Na [squares in Fig. 2(b)]. Cover¬ 
ing the Fe/CuN/Fe structures with single layer of Li- 
p(lxl) tends to reduce the value of lEC for N = 7 
and increases the absolute value of the lEC for antifero- 
magnetically coupled Fe/Cu 4 /Fe. Moreover, the mag¬ 
netic order is reversed from FM (Fe/Cus/Fe) to AF 
(Fe/Cus/Fe/Li) when the spacer thickness is five Cu lay¬ 
ers. Conversely, for Fe/Cug^s/Fe trilayers Li doping has 
the tendency of reducing the strength of AF coupling, 
and even the relative magnetic order in Fe/Cug/Fe is 
changed to FM coupling (for Fe/Cug/Fe/Li), as result 
of Li doping. In Fig. 2(b)] we have also presented the 
lEC for Fe/CuN/Fe/Li-c(2x2) multilayer, shown in cir¬ 
cles. Here reducing the percentage coverage of Li ions 
from p(lxl) to c(2x2) can only change the strength of the 
exchange coupling. Only for the case of Fe/Cuy/Fe/Li- 
p(lxl) the mutual magnetic order can be switch to AF 


coupling by reducing Li coverage [Li-c(2x2)] of the mag¬ 
netic surface. Interestingly, the change on the lEC of 
FeCu multilayers that are caused by p(lxl) configura¬ 
tion of both Na and Li ions are alike. Moreover, FeCu 
multilayers that are doped with Na-c(2x2) and Li-c(2x2) 
also show similar exchange coupling in most cases. All 
in all, except few cases the trends in which Li or Na 
ions change the lEC is analogous to one another, im¬ 
plying that similar tuning of the relative magnetization 
in FM/NM/FM multilayers can also be accomplished by 
other alkali metals. 


V. MAGNETIC HARDENING OF FePt 
MULTILAYERS WITH Na IONS 

We have demonstrated a remarkable effect of direct 
surface charging and alkali-metal-doping on the inter¬ 
layer exchange coupling. It is also worthwhile to inves¬ 
tigate the impact of Na ions on the magnetic anisotropy 
energy (MAE) which determines the stability of magnetic 
systems. In this regard materials with FePt composition 
are known to have high MAE*^’^^’"*® and are used for the 
current high density recording media^®. Here, we show 
that MAE of FePt multilayer can be further enhanced 
by adsorption of Na ions. In order to verify such way 
of magnetic hardening effect, we consider a model sys¬ 
tem which consists of two Fe layers separated with Pt 
multilayers and theses Fe/PtN/Fe stack is supported by 
Pt(OOl) substrate. We calculate the MAE as the differ¬ 
ence in energy between two axis of magnetization of the 
system, i.e. MAE = Ex — E^. Identical calculation 
is carried out for Fe/PtN bilayer on Pt(OOI), in order to 
to estimate the layer resolved MAE and its relation with 
the interlayer magnetic coupling®®. 

In Fig. 3 (a) the MAE of Fe/PtN/Fe trilayer is plotted 
as function of the Pt spacer thickness, with (squares) and 
without (circles) the adsorption of Na ions. The values 
of the MAE are presented for the stable magnetic config¬ 
urations, and in all cases the easy axis of magnetization 
is aligned perpendicular to the surface. The FePt mul¬ 
tilayers with smaller spacer thicknesses (N = 2,3) show 
high value of MAE relative to the thicker ones. A de¬ 
crease in the value of the MAE of Fe/PtN/Fe is observed 
when the Pt spacer thickness increases from N = 3 to 
N = 4. Further increment of the spacer thickness brings 
small change to the magnetic anisotropy and an aver¬ 
age MAE of 1 meV is observed. As shown in Fig. 3 
(b) (circles), for Fe/Pt 2 , 5 /Fe the stable magnetic order 
is found to be ferromagnetic, whereas antiferromagnetic 
coupling is observed for N = 3, 4 and 6. Indeed, the 
interlayer magnetic order for Fe/Pt 2 /Fe is in agreement 
with similar other systems®®. Additionally, experimental 
and theoretical studies have verified the existence of such 
variation of the interlayer exchange coupling between FM 
and AF coupling in similar systems®*’®^. 

Basically, Pt is very close to Stoner criterion and can 
be magnetized when it hybridizes with a ferromagnet. 
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FIG. 3. (Color online) (a) The magnetic anisotropy en¬ 
ergy (MAE) of the stable magnetic orders in Fe/PtN/Fe (cir¬ 
cle) and Na/Fe/PtN/Fe (square) for different thickness of Pt 
spacer is presented. The arrows in the square and circle signs 
indicate the respective magnetic order for which the magnetic 
anisotropy is evaluated, (b) The corresponding plot for the 
interlayer exchange coupling where positive and negative val¬ 
ues of the lEC represent ferromagnetic and antiferromagnetic 
couplings, respectively. 


For example, we found that the nearest Pt layer has a ^ 
magnetic moment up to 0.35 fiB, in Fe/PtN/Fe multi- 
layers. In such systems the competition of the induced gj 
magnetization from the upper and lower Fe layers deter¬ 
mines the magnetic configuration in the neighboring Pt 
layers, more importantly for N < 3, and probably the 
magnetic order within the multilayer. As a consequence 
of such competition the Fe/Pt 2 , 3 /Fe trilayers could be 
driven away from Stoner instability, leading to antifero- 
magnetic order. Certainly, the effect of the spin polarized 
QWS on the lEC still exits, specially for FePt multilayers 
with thicker Pt spacer. Even more such magnetic cou¬ 
pling can be affected by Na capping layer. For instance, 
strong change of the lEC is depicted for Fe/Pts^e/Fe, in 
Fig. 3 (b), by Na adsorption. Whereas Na capping layer 
induces moderate change for other FePt multilayers. Ad¬ 
ditional perspective on the tuning of the lEC in relation 
to the MAE and Na overlayer are pointed out in the latter 
part of Sec. VI. 

As already mentioned in the earlier sections one of the 
main goals of this study is to point out the impact of al¬ 
kali metal deposition on the magnetic anisotropy of FePt 
thin films. In Fig. 3 (a), contrary to the uncapped mul¬ 
tilayers, the FePt multilayers that are capped with Na 
monolayer show high MAE for thicker spacer thicknesses. 
Clearly one can see that Na adsorption strongly enhances 
the MAE when the Pt spacer thicknesses of four, five 
and six atomic layers. As an example, deposition of Na 


ions on Fe/Ptg/Fe trilayers (AF coupling) significantly 
enhances the MAE from ~ 1 meV to ^ 3.4 meV, in¬ 
creasing the energy barrier by more than three fold. Ex¬ 
ceptionally, for the system whereby the two Fe layers 
are interspaced with five Pt layer, Na ion capping in¬ 
crease the MAE to ~8 meV. Moreover, we have com¬ 
pared the change in the MAE and the lEC and it is 
clearly seen that the systems with the strongest lEC cor¬ 
responds to the highest MAE. Specifically, we have al¬ 
ready shown that Na deposition strongly increases the 
lEC of Fe/Pts/Fe trilayers from 1.5 meV to 22 meV. 
The next higher changes in the lEC induced by Na over¬ 
layer are observed when the spacer thicknesses of four and 
six Pt layers are employed. Likewise, the change in the 
MAE is relatively large. These results clearly hint that 
the enhancement of the MAE is almost directly related 
with the absolute value of the change in the interlayer ex¬ 
change coupling. Moreover, it verifies that the magnetic 
order and anisotropy of such complex alloy magnetic sys¬ 
tems, could be altered as result of ionic doping, leading 
to a change on the magnetic stability and dynamics^®’®®. 
It is experimentally proven fact that the magnetic or¬ 
der of such FM/NM/FM trilayers remains unchanged for 
thicker ferromagnetic layers®^’®^, and hence the interplay 
between the magnetic order and anisotropy. 
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FIG. 4. (Golor online) (a) The hysteresis loop of Fe/Pts/Fe 
trilayers with (full lines) and without (dotted lines) Na cap¬ 
ping layers. The simulations are performed for out-of-plane 
direction which is the easy axis of magnetization, (b) An illus¬ 
trative diagram of the exchange coupling parameters within 
the multilayer. 

Apart from the magnetic properties discussed above 
the magnetization reversal process is studied by simulat¬ 
ing the hysteresis loop for FePt multilayers. In order to 
to do this we apply an external magnetic field within the 
range of -12 to 12 T. The calculation is performed within 
the stochastic Landau-Lifshitz-Gilbert (LLG) equation. 
The LLG equation (Eq. 2) is integrated by using the 
Heun integration scheme®^, with a time step of At = 1 fs. 
The simulations temperature is set to be 5 K. Whereas 
the damping parameters A = 1.0 is used for the sake 
of computational efficiency. The values of the exchange 
coupling parameters, Ji, J 2 and J 3 , as well as the the es- 
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timated values of the MAE are presented in TABLE I®®. 
We have created model system of two Fe layers separated 

TABLE I. The exchange parameters and J 3 

(meV/(/rB)^)], spin magnetic moments Fei(Fe 2 )(;iB)] and 
the estimated MAE (in meV) of two Fe layers interspaced with 
five layers of Pt. 



Ji 

1 ■’^2 

J 3 

MAEi 

IMAE 2 

|m^ [Fei(Fe2)] 

Fe2/Pt5/Fei 

5.1 

4.9 

0.2 

0.1 

1 1-2 

1 12.9(2.7)1 

Fe2/Pt5/Fei/Na 

4.7 

4.9 

1.4 

1.6 

1 6.3 

1 12.7(2.9)1 


with an effective spacer. Fig. 4 illustrates the hystere¬ 
sis loops for both FePt multilayers with (full lines) and 
without (dotted lines) Na deposition. Both multilayers 
show square like hysteresis loops. One can see that the 
Na capping layer significantly increase the coercive field 
to Ha = ~8.5 T. Therefore, the atomistic spin model 
simulations clearly demonstrate that alkali metal doping 
of FePt multilayers favor domain nucleation at elevated 
temperature by increasing the energy barrier. Further¬ 
more it indicates that magnetic hardening effect can be 
accomplished by Na ion deposition. 

VI. ANALYSIS AND DISCUSSION OF THE 
RESULTS 

The findings that have been discussed until now are 
supported with analysis on the electronic structure. It 
is known that the interaction between the magnetic lay¬ 
ers is mediated by the confined delocalized electrons in 
Cu layers^®’®®, which is a direct consequence of the spin- 
polarized potential at the Fe/Cu interfaces. Indeed, it is 
observed that the net charge carriers for negative (posi¬ 
tive) charge doping spills (depletes) mainly at the surface, 
and minimum variation exists below the surface. At this 
stage there are two factors that can impose changes on 
the potential of the quantum-well like structure, where 
the first one is related to the changes in the boundary 
conditions as a result of the relative alignment of the two 
magnetic layers, parallel or anti-parallel. Secondly ad¬ 
dition or removal of extra charge can alter the Fe/Cu 
interface as well as the reflectivity of the confined states 
at the boundaries®®. These two modification will affect 
the nature of the confinement, being additional scat¬ 
tering constraint of the quantum-well states. In Fig. 5 
the effect of these two parameters on the spin polar¬ 
ized quantum-well states (QWSs) of Fe/Cu 4 /Fe trilayer 
is depicted. In these figures the s+p density of states 
of the coupling medium, for all Cu layers, are plotted 
for FM (a) and AF (b) coupling. We have considered 
three charge states, namely positive (<?+), neutral (go) 
and negative (g-). An electron or hole doping of the 
order of 1 e {h) is used for the charged systems. One 
can infer that the shift of the spin polarized QWS in the 
energy scale depends on the polarization of the charge 
doping and consequently this can alter the nature of the 


magnetic coupling. For both magnetic orders the pos¬ 
itive (negative) charge doping shifts the QWS towards 
higher (lower) energies, on the other hand, the variation 
of the magnetic order between the magnetic layers can be 
directly inferred from the change in the integrated den¬ 
sity of state (DOS) which significantly contributes to the 
total energy®®’®"* , E = n(s)(e — sp)ds, where n(e) is 
the density of states. Thus, the change in the lEC is ex¬ 
plained from the shift in the DOS of the magnetic layer 
or the coupling medium for FM and AF configuration at 
different charge states. As shown in Fig. 5 the shift of 
the QWS induced by positive charging is more enhanced 
for AF coupling. The dependence of such band shift on 
the magnetic order can lead to stable ferromagnetic cou¬ 
pling in Fe/Cu 4 /Fe for positive charging, as it is found 
in section III. It has to be noted that the QWS are cal¬ 
culated by sampling the k-points near the center of the 
Brillouin zone. In order to get additional information we 
also analyzed the total density of states. 



FIG. 5. The spin polarized quantum well states of all Cu 
spacer layers for FM (a and c) and AF (b and d) orders, 
at different charge states of the supercell are depicted. The 
majority and minority DOS are presented on the upper (a 
and b) and lower (c and d) panel, respectively. The density of 
states are plotted for positive (g+=l h, blue full line), negative 
(g_=l e, red dashed line) and neutral (dotted line) Fe/Cu 4 /Fe 
multilayers, close to the center of Brillouin zone. The vertical 
dashed line at zero energy level represents the Fermi energy. 

Studying the variation of the electronic structure of 
the top FeCu interface can also provide more insight in 
the change of the lEC caused by surface charging. First, 
we would like to focus on the electronic structure of the 
magnetic layer in order to elucidate the magnetic rever¬ 
sal scenario once more. For both magnetic orders the 
shift of the minority bands, caused by the surface charg¬ 
ing process, show minimum difference and are omitted. 
Unlike for the case of Cu-QWS (Fig. 5) the majority and 
minority bands of the Fe layer show opposite shift rela¬ 
tive to the neutral system. This leads to an increase or 
decrease in the exchange splitting as well as the magnetic 
moment. As shown Fig. 6(a) and (b) the majority QWS 
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of Cu spacer is shifted to higher energy as Fe/Cu 4 /Fe tri¬ 
layer is positively charged but the majority DOS of Fe are 
shifted to lower energies. This is caused by the the net 
charge accumulated/depleted and the reflectivity of the 
confined electrons within the quantum-well like structure 
(Fe-Cu-Fe) can be modified therewith. 

In Fig. 6 (a) and (b) the majority DOS for d-orbitals of 
the top Fe layer, for parallel and anti-parallel alignment, 
at two different charge states is presented. Similarly, the 
density of states of the Cu atoms that are interfaced with 
Fe layer is plotted in Fig. 6 (c) and (d). From these 
density of states one can infer that there is a hybridiza¬ 
tion between the sp states of Cu with d-orbitals of Fe. 
Explicitly, injecting positive charge of 1 /i per unit cell 
induces a strong shift (cr) of the majority d-bands of Fe 
towards the lower energies by cr = 140 meV and cr = 65 
meV for FM and AF configurations, respectively. This 
implies that the majority bands of Fe, in the case of fer¬ 
romagnetic coupling, is more sensitive to positive charge 
doping which leads to a lower total energy and stable 
ferromagnetic coupling. The reason for such change of 
the majority bands is related to the sp — d hybridiza¬ 
tion of the Cu layer with top Fe layer. When the system 
is either negatively or positively charged the sp states 
are affected substantially and it leads to a change in the 
sp — d hybridization. Evidently the corresponding shift 
of the DOS of Cu induced by positive charge doping, to 
the lower energies, has similar tendency like that of the 
majority bands of Fe, see Fig. 6 (c) and (d). 



FIG. 6. The majority electron density of states (for d-orbitals) 
of the top Fe layer Fe/Cu 4 /Fe, for both FM (a) and AF (b) 
configurations are shown. The shift in the DOS (cr) are com¬ 
pared for the neutral system and when it is charged with 1.0 h 
per unit cell. The majority density of states for the Cu layer 
interfaced with top Fe layer is shown in the lower panel. The 
Cu-DOS is plotted for the cases where Fe/Cu 4 /Fe trilayer is 
coupled ferromagnetically (c) and antiferomagnetically (b). 
The density of states are integrated over the whole Brillouin 
zone. The vertical dashed line at zero energy level represents 
the Fermi energy. 

Similar to the case of direct surface charging some in¬ 
sights on the variations of the magnetic order (by alkali 
metal doping) can be inferred by comparing the QWS 
and the integrated density of states. However, the mech¬ 
anism in which the lEC is changed by Na or Li doping in¬ 


volves both chemical and electronic effects, unlike that of 
direct surface charging. In Fig. 7 the minority quantum 
well states within the coupling medium of Fe/Cu 4 /Fe are 
calculated near the center of Brillouin zone and presented 
for FM (a) and AF (b) coupling. In this figure the QWS 
for both FM and AF orders are displaced to the lower 
energy as result of Na doping, which is analogous to the 
case of negative surface charging. Additionally, in Fig. 7 
the minority DOS of the top Fe layer in Fe/Cu 4 /Fe tri¬ 
layer is also shown for both FM (c) and AF (d) couplings 
as well as the DOS of Na capping layers. In both plots, 
i.e. upper panel (Cu-QWS) and lower panel (Fe-DOS) 
the minority bands are correspondingly shifted to the 
lower energy as result of Na capping layer. Clearly the 
effect is more pronounced for the case of FM coupling 
than the AF ones. Following Na adsorption more states 
are displaced to the lower energies in the case of FM cou¬ 
pling than AF one, leading to stable FM order. Indeed, 
similar changes are also observed for the bottom mag¬ 
netic layer which is mediated by quantum well states. In 
section IV, we have seen that either reducing the percent¬ 
age coverage of Na [Na-c(2x2)] or Li-p(Ixl) can reverse 
the magnetic order of Fe/Cu 4 /Fe from AF order to FM 
order. To get some insight about these differences we 
have added the minority-DOS plot for top Fe layer [grey 
dashed line, Fig. 7 (c) and (d)] when the Fe/Cu 4 /Fe is 
doped with Li-p(lxl). It seems that Li does not bring 
strong change to the spin polarized electronic structure 
of the Fe layer which can induce magnetic reversal. Ac¬ 
cording to Fig. 7 (c) and (d) the shift of the minority 
bands of Fe to the lower energies by Li-p(lxl) overlayer, 
for FM and AF couplings, is not as pronounced as in the 
case of Na-p(lxl) doping. In the contrary, more minor¬ 
ity DOS of Fe layers shifts to lower energy for the case 
of AF coupling which might only change the value of the 
exchange coupling. 

Once we have discussed the impact of varying the sur¬ 
face electronic structure with surface charging or alkali 
metal doping. Now, it is important to compare the effects 
of both surface charging techniques on the lEC. From 
Fig. 1 and Fig. 2 one can see that covering the magnetic 
layers with Na-c(lxl) does not necessarily lead to sim¬ 
ilar change like that of electron doping (1 e/unit cell). 
As an example, for Fe/Cu 4 /Fe trilayer electron doping 
and Na [Na-c(lxl)] deposition have opposite effect rela¬ 
tive to the neutral system. However, the impact of the 
submonolayer Na-c(2x2) and the direct surface charging 
techniques consistently favors the same magnetic orders 
for almost all Cu spacer thickness. 

Using atomistic spin model simulations, perhaps for 
the first time we revealed magnetic hardening mechanism 
in FePt multilayers by Na ion doping. These variations in 
the MAE are mainly caused by a combined effect of the 
interlayer exchange coupling and the changes in surface 
electronic structure induced by Na. In order to eluci¬ 
date the first case we compared the MAE for three dif¬ 
ferent cases, i.e. Fe/Pts, Fe/Pts/Fe and Fe/Pts/Fe/Na 
on Pt(OOl). The MAE of the first structure is found to 


















Energy (eV) 


FIG. 7. The effect of Na capping layer in the confined states 
of the quantum well like structure, Fe/Cu 4 /Fe. The DOS 
of the QWS are presented for minority states with [(a), full 
line] and without Na doping [(b), dashed line]. The minor¬ 
ity electronic density of states (for d-orbitals) of Fe depicted 
without (full lines) and with (dashed lines) adsorption of Na 
ions. All the DOSs are plotted for the top Fe layer within 
the FM (c) and AF (d) configurations of Fe/Cu 4 /Fe. The 
dashed gray line shows the minority bands of the top Fe in 
Fe/Cu 4 /Fe/Li-p(lxl) multilayer. The DOS for the Na layer 
is also shown in dotted lines and it is multiplied by factor of 
ten. Since the shift of the majority DOS as result of Na ion 
doping is relatively similar, for both FM and AF order, and 
are not included. The vertical dashed line at zero energy level 
represents the Fermi energy. 


be 1.6 meV. When the bilayer Fe/Pts is terminated with 
Fe layer (Fe/Pts/Fe), h leads to the onset of the effect of 
interlayer exchange coupling on the magnetic anisotropy. 
The added (top) magnetic layer which couples ferromag- 
netically with the bottom Fe layer brings small change 
to the MAE. However, deposition of Na ion increases the 
MAE enormously. Here the strong change of the MAE 
is the ramifications of the lEC and the change in the 
electronic structure of the FePt interfaces. This implies 
that both the magnetic order and the interface electronic 
structure are relevant parameters that could possibly af¬ 
fect the MAE by varying the spin-orbit coupling. Es¬ 
sentially, the contribution of the spin-orbit coupling to 
each axis of magnetization, in-plane or out-of plane, de¬ 
termines the MAE. For instance, in FePt multilayers the 
coupling between occupied and unoccupied states of ds 
and d,r orbitals are identified as the dominant ones among 
the minority d—orbitals'^^’®'^, where ds = dx 2 -y 2 + dxy 
and = dxz{yz) + For the case of Fe/Pts/Fe/Na 
the orbitals ds (d^) are less (more) degenerate near to 
the Fermi energy for the top Fe layer as compared to 
the lower ones. This tends to lower the energy when the 
magnetization of the multilayer is perpendicular to the 
surface, leading to a substantial increase in the MAE for 
the bottom FePt interfaces. In addition the magnetic 
moment as well as the induced magnetic moment for the 
lower FePt interfaces has increased as result of Na cap¬ 
ping, and conversely reduces for top layers. 

On the other hand some of the FePt multilayers with 
high MAE have stable antiferomagnetic coupling, imply¬ 
ing that these FePt thin films could be essential can¬ 


didates for AF domains®^. Besides, AF domains have 
highly minimized dipolar magnetic interaction between 
the neighboring bits, and could play a crucial role in 
reducing the density of the current magnetic storage 
device"^®. Additionally such AF multilayers can be an 
integral part of an exchange biased system and could 
have significant impact in the quest of ultra-fast mag¬ 
netic switching®®’®®. In fact, the switching procedure of 
the relative magnetization by incorporating alkali-metal 
doping can also be perceived as follows. Regardless of the 
original magnetic order in FM/NM/FM trilayer, first one 
has to overcome the energy barrier (MAE) of the top fer- 
romagnet in order to switch to the opposite direction. 
The magnetic anisotropy of each magnetic layer in FePt 
multilayer is estimated in the earlier section. This en¬ 
ergy barrier can certainly be altered by introducing alkali 
metal in to the system, which varies the interface elec¬ 
tronic structure. The charge rearrangements are then 
transformed to magnetic anisotropy through spin-orbit 
coupling which results in changes of the relative orienta¬ 
tion of the atomic moments. But both surface charging 
techniques does not always lead to magnetization rever¬ 
sal. Instead, a favorable conditions will be created for 
other means of spin manipulation mechanisms, e.g. spin 
transfer torque^’® or electrostatic gating®^’®®’®®, leading 
to switching mechanism assisted with alkali metal dop¬ 
ing. Therefore, once the energy barrier is completely 
overcome the system probably settles down to the re¬ 
versed state. 


VII. CONCLUSION 

In summary, we have presented a comprehensive study 
which reveals the effect of two surface-interface charg¬ 
ing techniques on the lEC of magnetic layers across 
Cu and Pt spacers as well as their MAE. Explicitly, in 
some cases the mutual magnetic orders of fully metal¬ 
lic FM/NM/FM-like structure can be remarkably tuned 
between ferromagnetic and antiferromagnetic states by 
addition or removal of extra charge. Moreover, employ¬ 
ing alkali metals as dopants we have shown a natural way 
of space charge alteration mechanism which can reverse 
the relative magnetic orders in alloyed multilayers and 
strongly enhance the MAE of Fe/Pt. These effect leads 
us to the conclusion that alkali-metal-induced charging 
magnetic surfaces can assist the process of magnetiza¬ 
tion reversal and also controlling the spin states as well 
as the spin dynamics. All in all, our findings suggest 
that involving alkali metals in field of spintronics could 
have significant impact in order to control and design new 
functional spintronic devices or quantum computing. 
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